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ABSTRACT 

We report an analysis of new V and I CCD time-series photometry of the distant globular 
cluster NGC 6229. The principal aims were to explore the field of the cluster in search of 
new variables, and to Fourier decompose the RR Lyrae light curves in pursuit of physical 
parameters.We found 25 new variables: 10 RRab, 5 RRc, 6 SR, 1 CW, 1 SX Phe, and two 
that we were unable to classify. Secular period changes were detected and measured in some 
favourable cases. The classifications of some of the known variables were rectified. 

The Fourier decomposition of RRab and RRc light curves was used to independently 
estimate the mean cluster value of [Fe/H] and distance. From the RRab stars we found 
[Fe/H]uvES=—1-31 ± 0.01 (statistical) ± 0.12(systematic) ([Fe/H]zw = —1.42),and a dis¬ 
tance of 30.0 ± 1.5 kpc, and from the RRc stars we found [Fe/H] uves=— 1-29 ±0.12 and 
a distance of 30.7 ± 1.1 kpc, respectively. Absolute magnitudes, radii and masses are also 
reported for individual RR Lyrae stars. Also discussed are the independent estimates of the 
cluster distance from the tip of the RGB, 34.9±2.4 kpc and from the P-L relation of SX 
Phe stars, 28.9±2.2 kpc. The distribution of RR Lyrae stars in the horizontal branch shows a 
clear empirical border between stable fundamental and first overtone pulsators which has been 
noted in several other clusters; we interpret it as the red edge of the first overtone instability 
strip. 

Key words: globular clusters: individual (NGC 6229) - stars:variables: RR Lyrae, SX Phe, 
SR 


1 INTRODUCTION 

As part of a long-term programme of CCD time-series photom¬ 
etry of globular clusters (GC) aimed at updating their variable 
star census and offering independent estimates of their metallic- 
ity and distance, in the present paper we carry out the study of 
the globular cluster NGC 6229 (C1645+476 in the lAU nomencla- 
ture) (a = 16‘’46”58.6®, <5 = +47°31'36.4", J2000; I = 73.64°, 
b = +40.31°). This is a distant outer-halo cluster (d ~ 30 kpc) 
which is moderately concentrated and hence variable stars have 
been identified in regions very near its core since very early in pho¬ 
tographic studies of the cluster. The Catalogue of Variable Stars 
in Globular Clusters (CVSGC; Clement et al. 2001; 2013 edition) 

* Based on observations collected with the 2.0-ni telescope at the Indian 
Astronomical Observatory. 


lists 48 variables, mostly of the RR Lyrae type (RRL). The first 
variable star in NGC 6229, a type II cepheid, was announced by 
Davis (1917) and later named as V8 by Baade (1945) who, us¬ 
ing 46 photographic plates obtained at Mount Wilson in the years 
1932-1935, discovered another twenty variables, V1-V7, V9-V21. 
Sawyer (1953) added the W Virginis or CW star V22 indicating that 
its likely period is longer than one day. Nearly forty years passed 
before two more variables were found by Carney et al. (1991), al¬ 
ready with a CCD detector, and later confirmed by Borissova et 
al. (2001) (BCVOl); these are the two RRc stars V23 and V24. 
Spassova & Borissova (1996) suggested 12 possible variables, only 
seven of which were confirmed by BCVOl; nevertheless, these au¬ 
thors numbered all 12 stars as V25-36. However, the variability 
status and/or type of stars V25, V26, V27, V28, V29 and V30 have 
remained dubious in the CVSGC. In their paper, BCVOl also re¬ 
ported 12 new variables, V37-48. 
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In most contemporary variable star studies of NGC 6229 using 
CCD photometry, the recipes based on PSF or aperture photome¬ 
try (Stetson 1987) and image subtraction algorithms such as ISIS 
(Alard 2000) have been used. In the present paper we employ the 
DanDIA implementation of difference image analysis (DIA) to ex¬ 
tract high-precision photometry for all of the point sources in the 
field of NGC 6229. In the 13 years elapsed since the last CCD 
time-series analysis of the cluster (BCVOl), important improve¬ 
ments have taken place in the concept of DIA, which are described 
in detail by Bramich (2008) and Bramich et al. (2013). Hence, we 
considered that a new time-series study may reveal a useful update 
of the variable star population in NGC 6229. 

In previous papers we have demonstrated how DIA performs 
exceptionally well in finding new variables and low signal-to-noise 
effects in light curves such as amplitude modulations in RR Lyrae 
stars, even in the crowded central regions of globular clusters (e.g. 
Lazaro et al. 2006; Bramich et al. 2011; Arellano Ferro et al. 2011, 
2013 and references therein). In the present paper we report the re¬ 
sults of the analysis of new time-series photometry in the V and I 
filters of NGC 6229. We have Fourier decomposed the RR Lyrae 
(RRL) light curves to calculate [Fe/H], distance, effective tempera¬ 
ture, mass and radius for individual stars. Mean [Fe/H] and distance 
are hence good estimates of the parent globular cluster values. 

The plan of the paper is as follows. In § 2 we describe the 
observations and data reduction. In § 3 the approaches to new vari¬ 
able star identification are described and we report the new vari¬ 
ables found and discuss the secular period changes calculated in 
some favourable cases. § 4 contains the Fourier light curve decom¬ 
position of some of the RR Lyrae stars and we report their cor¬ 
responding individual physical parameters. § 5 is dedicated to the 
discussion of the distribution of RRL stars in the HB in NGC 6229 
and in other recently studied clusters. In § 6 the metallicity and 
distance of the parent cluster are inferred from the RRL stars. In¬ 
dependent distance calculations are also carried out from the P-L 
relation of SX Phe stars, and the luminosity of the tip of the red 
giant branch (RGB). In § 7 we summarise our conclusions. Finally, 
the classification and peculiarities of individual stars are discussed 
in Appendix A. 


2 OBSERVATIONS AND REDUCTIONS 
2.1 Observations 

The observations were performed on 20 nights between 7th April, 
2010 and 19th June, 2013 with the 2.0 m telescope at the Indian 
Astronomical Observatory (lAO), Hanle, India, located at 4500 m 
above sea level. A total of 346 and 354 images were obtained in the 
Johnson-Kron-Cousins V and / filters, respectively. The detector 
was a Thompson CCD of 2048x2048 pixels with a scale of 0.296 
arcsec/pix, translating to a field of view (FoV) of approximately 
10.1x10.1 arcmin^. 

The log of observations is given in Table ^ where the dates, 
number of frames, exposure times and average nightly seeing are 
recorded. 


2.2 Difference Image Analysis 

We employed the technique of difference image analysis (DIA) to 
extract high-precision photometry for all of the point sources in the 


Table 1. The distribution of observations of NGC 6229 for each filter. 
Columns Ny and A/ give the number of images taken with the V and / 
filters respectively. Columns ty and tj provide the exposure time, or range 
of exposure times, employed during each night for each filter. The average 
seeing is listed in the last column. 


Date 

Ny 

ty (s) 

Nj 

ti (s) 

Avg seeing (”) 

20100407 

14 

500 

13 

100 

2.2 

20100505 

2 

500 

3 

100 

1.9 

20110411 

6 

500-600 

8 

100-200 

2.7 

20110412 

6 

400-500 

8 

100-200 

2.1 

20110413 

11 

350-500 

12 

75-150 

1.8 

20110414 

15 

350-400 

16 

70-80 

1.9 

20110609 

2 

600 

1 

200 

1.7 

20110611 

16 

450 

15 

120-200 

1.8 

20110612 

25 

110-450 

24 

110-120 

1.7 

20110709 

13 

140-600 

12 

140-160 

2.4 

20110710 

14 

400-600 

15 

80-120 

2.1 

20120428 

1 

600 

1 

200 

2.6 

20120429 

5 

600 

1 

200 

2.1 

20130304 

54 

30-220 

54 

12-80 

2.5 

20130305 

45 

40-100 

47 

15-45 

2.5 

20130306 

41 

50-70 

40 

18-30 

2.3 

20130416 

25 

240-400 

28 

30-150 

2.5 

20130417 

19 

400 

20 

125-150 

2.3 

20130618 

18 

300 

20 

80-100 

1.7 

20130619 

14 

300 

16 

100 

1.7 

Total: 

346 


354 
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Figure 1. Transformation relations in the V and I band-passes between the 
instrumental and the standard photometric systems using a set of standard 
stars in the field of NGC 6229 kindly provided by Peter Stetson. 

images of NGC 6229 and we used the DanD I pipeline for the 
data reduction process (Bramich et al. 2013). DanDIA includes an 
algorithm that models the convolution kernel matching the PSF of 
a pair of images of the same field as a discrete pixel array (Bramich 
2008). The DanDIA pipeline performs standard bias level and flat 
field corrections of the raw images, and creates a reference image 
for each filter by stacking a set of registered best-seeing calibrated 
images. 

^ DanDIA is built from the DanIDL library of IDL routines available at 
http://www.danidl.co.uk 
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Figure 2. The rms magnitude deviations as a function of the mean mag¬ 
nitudes V and I. Blue and green circles are used for RRab and RRc stars 
respectively. Other symbols are: red circle, V36 reclassified as RRc in this 
paper; orange circles, CWA, CWB or CW?; red triangles, SR; turquoise tri¬ 
angles, variable blue stragglers; red crosses, previously suspected variables 
not confirmed in this paper. 


We constructed two reference images, one in the V filter and 
another in I. For each of these reference images, 5 calibrated im¬ 
ages were stacked with total exposure times of 1800 and 370 sec¬ 
onds in the V and I filters, respectively, and with PSF FWHMs of 
~5.2 and ~4.0 pix, respectively. In each reference image, we mea¬ 
sured the fiuxes (referred to as reference fluxes) and positions of 
all PSF-like objects (stars) by extracting a spatially variable (with a 
third-degree polynomial) empirical PSF from the image and fitting 
this PSF to each detected object. The detected stars in each image 
in the time-series were matched with those detected in the corre¬ 
sponding reference image, and a linear transformation was derived 
which was used to register each image with the reference image. 

For each filter, a sequence of difference images was created 
by subtracting the relevant reference image, convolved with an ap¬ 
propriate spatially variable kernel, from each registered image. The 
spatially variable convolution kernel for each registered image was 
determined using bilinear interpolation of a set of kernels that were 
derived for a uniform 6x6 grid of subregions across the image. 

The differential fluxes for each star detected in the reference 
image were measured on each difference image. Light curves for 
each star were constructed by calculating the total flux ftot{t) in 
ADU/s at each epoch t from: 


— /ref + 


/diff(t) 

p{t) 


( 1 ) 


Figure 3. The top panel shows the colour-magnitude diagram of NGC 6229. 
Symbols and colours are as in Fig.[^ All variables labelled with numbers 
larger than 48 are new variables announced in the present work. See the 
text for a discussion on individual stars. The three dashed boxes are re¬ 
gions where individual star light curves were explored for variability. The 
isochrones are from the Victoria-Regina stellar models (VandenBerg et al. 
2014) for an age of 12 Gyr and [Fe/H]=—1.42 (red) and —1.31 (blue) re¬ 
spectively. The bottom panel is a blow-up of the HB region, where triangles 
represent stars with amplitude modulations and the two vertical black lines 
are two empirical loci of the border between RRab and RRc distributions 
seen in NGC 5024 (solid) and NGC 4590 (dashed) after considering the 
corresponding reddenings (see the discussion in § 0. Two ZAHB mod- 
els calculated from the Victoria-Regina stellar models for [Fe/H]=—1.42, 
[a/Fe]=+0.4 and two values of helium content Y 0.25 and 0.27 are included 
(see the discussion in § 5). 


where /ref is the reference flux (ADU/s), /diff(t) is the differential 
flux (ADU/s), and p{t) is the photometric scale factor (the integral 
of the kernel solution). Conversion to instrumental magnitudes was 
achieved using: 


mins(t) = 25.0 - 2.5 log [/tot(/)] , (2) 

where mins(t) is the instrumental magnitude of the star at time t. 
Uncertainties were propagated in the correct analytical fashion. 

The above procedure and its caveats have been described in 
detail in Bramich et al. (2011), to which the interested reader is 
referred for further details. 
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Table 2. Time-series V and I photometry for all the conhrmed variables in our held of view. Filter and epoch of mid-exposure are listed in columns 2 and 
3, respectively. The standard Mstd and instrumental mins magnitudes are listed in columns 4 and 5, respectively, corresponding to the variable star in column 
1. The uncertainty on mins is listed in column 6, which also corresponds to the uncertainty on Mstd- For completeness, we also list the quantities /ref, fdiff 
and p from Equation[T]in columns 7, 9 and 11, along with the uncertainties (Jref and in columns 8 and 10. This is an extract from the full table, which is 
available with the electronic version of the article (see Supporting Information). 


Variable 

Star ID 

Filter 

HJD 

(d) 

Mstd 

(mag) 

ru-ins 

(mag) 

(Tm 

(mag) 

fret 

(ADUs-1) 

(T’ref 

(ADUs-1) 

/diff 

(ADUs-1) 

(T’diff 

(ADUs-1) 

P 

VI 

V 

2455294.32981 

17.971 

19.087 

0.009 

179.343 

0.957 

+26.504 

0.931 

0.5048 

VI 

V 

2455294.34033 

18.017 

19.133 

0.007 

179.343 

0.957 

+26.841 

0.923 

0.6244 

VI 

I 

2455294.31704 

17.435 

18.422 

0.029 

328.424 

2.483 

+44.876 

4.482 

0.4514 

VI 

I 

2455294.32313 

17.432 

18.418 

0.023 

328.424 

2.483 

+52.696 

4.111 

0.5226 

V2 

V 

2455294.32981 

18.140 

19.251 

0.011 

194.212 

0.956 

+2.550 

0.944 

0.5048 

V2 

V 

2455294.34033 

18.147 

19.259 

0.009 

194.212 

0.956 

+2.271 

0.942 

0.6244 

V2 

I 

2455294.31704 

17.530 

18.514 

0.032 

336.249 

3.188 

+25.571 

4.560 

0.4514 

V2 

I 

2455294.32313 

17.552 

18.536 

0.026 

336.249 

3.188 

+25.545 

4.237 

0.5226 


2.3 Photometric Calibrations 

2.3.1 Relative calibration 

All photometric data suffer from systematic errors to some level 
that sometimes may be severe enough to be mistaken for bona fide 
variability in light curves. However, multiple observations of a set 
of objects at different epochs, such as time-series photometry, may 
be used to investigate, and possibly correct, these systematic er¬ 
rors (see for example Honeycutt 1992). This process is a relative 
self-calibration of the photometry, which is being performed as a 
standard post-processing step for large-scale surveys (e.g. Padman- 
abhan et al. 2008; Regnault et al. 2009). 

We apply the methodology developed in Bramich & Freudling 
(2012) to solve for the magnitude offsets Zk that should be applied 
to each photometric measurement from the image k. In terms of 
DIA, this translates into a correction (to first order) for the system¬ 
atic error introduced into the photometry from an image due to an 
error in the fitted value of the photometric scale factor p. We found 
that, for either filter, the magnitude offsets that we derive are of the 
order of ~ 0.04 mag with ~ 10% of worse cases reaching ~ 0.2 
mag. Applying these magnitude offsets to our DIA photometry im¬ 
proves the light curve quality, especially for the brighter stars. 


2.3.2 Absolute calibration 

Standard stars in the field of NGC 6229 were not included in the 
online collection of Stetson (2000 J^at the time of preparing this pa¬ 
per. However, Prof. Stetson kindly provided us with a set of prelim¬ 
inary standard stars which we have used to transform instrumental 
vi magnitudes into the standard V7 system. 

The standard minus the instrumental magnitudes show mild 
dependences on the colour, as can be seen in Figj^ The transfor¬ 
mations are of the form: 


Vstd = ^ + 0.112(±0.026)(u -i)- 1.177(±0.033), (3) 

7,,^ = i + 0.066(±0.018)(u - z) - 1.022(±0.023). (4) 

All of our VI photometry for the variable stars in the FoV of 
our collection of images of NGC 6229 is provided in Table Just 
a small portion of this table is given in the printed version of this 
paper, while the full table is only available in electronic form. 

Fig.j^shows the rms magnitude deviation in our V and I light 
curves, after the application of the relative photometric calibration 
of Section [?.3.1| as a function of the mean magnitude. 

To help us discuss the variable star search and classifica¬ 
tions, we have built the colour-magnitude diagram (CMD) of Fig. 
I^by calculating the inverse-variance weighted mean magnitudes 
of 2113 stars with V and I magnitudes. For better precision, the 
periodic variables like RRL and SX Phe are plotted using their 
intensity-weighted mean magnitudes < V > and colours < V > 
— < I >. The use of < y — 7 > instead does not produce a sig¬ 
nificant difference in the plotting on the CMD as they differ from 
< V > — < I > within 0.01 mag. The regions marked on the 
plot, individual stars, and the HB inset will be discussed below in 
the corresponding sections. 


2.4 Astrometry 

A linear astrometric solution was derived for the V filter reference 
image by matching 84 hand-picked stars with the UCAC4 star cata¬ 
logue (Zacharias et al. 2013) using a field overlay in the image dis¬ 
play tool GAI (Draper 2000). We achieved a radial RMS scatter 
in the residuals of ~0.2 arcsec. The astrometric fit was then used to 
calculate the J2000.0 celestial coordinates for all of the confirmed 
variables in our field of view (see Table [^. The coordinates corre¬ 
spond to the epoch of the V reference image which pertains to the 


^ http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/coriimunity/STETSON/standards ^ http://star-www.dur.ac.uk/~pdraper/gaia/gaia.html 
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Table 3. General data for all of the confirmed variables in NGC 6229 in the FoV of our images. Stars V49-V73 are new discoveries in this work. Label ’BF is 
used for Blazhko variables. The number of digits quoted for < y > and < / > depend on the scatter and/or amplitude modulations seen in each light curve. 


Variable 

Star ID 

Variable 

Type 

< y > 

(mag) 

< / > 
(mag) 

Av- 

(mag) 

Ai 

(mag) 

P (days) 
this work 

/3 

(d Myr-i) 

HJDmax 
(+2 450 000) 

P (days) 
(BCVOl) 

RA 

(J2000.0) 

Dec. 

(J2000.0) 

VI 

RRab 

17.936 

17.450 

1.16 

0.82 

0.585700 

+0.120 

6399.3193 

0.5857 

16:46:56.34 

-^47:29:53.3 

V2 

RRab-Bl 

18.07 

17.56 

1.01 

0.63 

0.555239 

+0.044 

6357.4046 

0.5553 

16:46:51.62 

-^47:31:43.5 

V3 

RRab 

18.13 

17.59 

0.89 

0.67 

0.575259 

-0.432 

5753.1692 

0.5752 

16:46:39.38 

-^47:32:19.3 

V4 

RRab-Bl 

18.09 

17.53 

0.86 

0.58 

0.566196 

-0.500 

5664.3882 

0.5662 

16:46:53.12 

+47:31:24.3 

V5 

RRab 

18.160 

17.605 

1.05 

0.68 

0.533602 

-0.056 

5725.3179 

0.5336 

16:47:00.11 

+47:32:23.0 

V6 

RRab-Bl 

18.08 

17.54 

0.85 

0.52 

0.559349 

- 

6462.3056 

0.5593 

16:47:03.03 

+47:32:20.7 

V7 

RRab-Bl 

18.02 

17.55 

1.15 

0.77 

0.506969 

- 

5294.3403 

0.5016 

16:46:54.60 

+47:30:48.8 

V8 

CWA 

15.50 

14.62 

0.75 

0.58 

14.846 

- 

6462.3310 

14.8405 

16:46:58.33 

+47:30:56.9 

V9 

RRab-Bl 

17.95 

17.46 

0.86 

0.63 

0.542841 

-0.204 

6357.4281 

0.5428 

16:46:54.85 

+47:32:17.0 

VIO 

RRab 

18.101 

17.533 

1.00 

0.65 

0.554738 

-0.280 

6463.3279 

0.5547 

16:46:55.76 

+47:32:51.4 

Vll 

RRab 

18.119 

17.566 

0.92 

0.59 

0.563226 

- 

5752.2791 

0.5632 

16:47:01.08 

+47:31:14.1 

V12 

RRab-Bl 

17.91 

17.67 

1.16 

0.95 

0.490812 

- 

6463.3125 

0.4908 

16:47:02.11 

+47:31:15.4 

V13 

RRab 

18.106 

17.591 

0.93 

0.62 

0.547345 

+0.056 

6357.3728 

0.5473 

16:47:12.51 

+47:32:41.0 

V14 

RRab 

17.840 

17.353 

1.04 

0.66 

0.465926 

+0.096 

5725.3366 

0.4659 

16:46:57.21 

+47:30:48.1 

V15 

RRc 

18.12 

17.72 

0.52 

0.32 

0.271379 

+0.024 

5722.2369 

0.2714 

16:47:02.07 

+47:32:06.7 

V16 

RRc-Bl 

17.98 

17.61 

0.49 

0.30 

0.322836 

-fO.780 

6399.2877 

0.3227 

16:47:03.36 

+47:31:15.0 

V17 

RRc-Bl 

17.91 

17.59 

0.46 

0.35 

0.325141 

- 

6463.2279 

0.3252 

16:46:49.25 

+47:30:23.4 

V18 

RRab 

18.066 

17.497 

0.79 

0.51 

0.579039 

- 

5725.3594 

0.5791 

16:46:55.13 

+47:32:11.0 

V19 

RRab 

18.130 

17.631 

1.29 

0.87 

0.475960 

-0.024 

6462.3457 

0.4759 

16:47:03.99 

+47:30:54.9 

V20 

RRab-Bl 

18.17 

17.68 

1.30 

0.88 

0.465952 

-0.232 

6046.3548 

0.4660 

16:46:56.03 

+47:31:02.8 

V21 

RRab 

18.063 

17.546 

0.98 

0.66 

0.564426 

- 

6462.3941 

0.5644 

16:47:10.32 

+47:30:38.0 

V22 

CWA 

15.28 

14.43 

0.88 

0.72 

15.846 

- 

6047.4139 

15.8373 

16:46:58.95 

+47:31:28.1 

V23 

RRc 

17.868 

17.433 

0.38 

0.22 

0.397688 

- 

5722.2369 

0.3977 

16:46:56.72 

+47:32:36.0 

V24 

RRc 

18.016 

17.668 

0.44 

0.30 

0.303441 

- 

5752.2182 

0.3034 

16:46:57.45 

+47:30:40.5 

V27 

CW? 

16.83 

16.35 

0.49 

0.31 

- 

- 

6358.3593 

1.13827 

16:46:59.85 

+47:31:46.1 

V29 

RRab-Bl 

18.30 

17.44 

1.50 

0.74 

0.629509 

- 

6357.4207 

- 

16:47:00.39 

+47:31:41.5 

V30 

RRab 

18.7 

17.9 

0.70 

0.35 

0.498700 

- 

6356.3490 

- 

16:46:58.63 

+47:31:23.3 

V31 

RRab 

17.77 

17.24 

0.88 

0.61 

0.537851 

- 

6462.3941 

0.6989 

16:46:59.30 

+47:31:18.9 

V32 

RRab-Bl 

18.17 

17.63 

0.61 

0.49 

0.603805 

- 

5752.2266 

0.3765 

16:46:58.34 

+47:31:17.7 

V33 

RRab 

18.26 

17.60 

1.24 

0.78 

0.517500 

- 

6356.3640 

0.5176 

16:46:58.24 

+47:32:01.1 

V34 

RRab 

18.14 

17.58 

0.99 

0.61 

0.555409 

- 

6356.3936 

0.5554 

16:47:01.43 

+47:31:44.7 

V35 

RRab 

18.062 

17.501 

0.31 

0.23 

0.634460 

- 

5665.3884 

0.6345 

16:47:01.53 

+47:31:57.2 

V36 

RRc 

17.942 

17.569 

0.32 

0.22 

0.264133 

- 

6399.2877 

0.2641 

16:46:56.51 

+47:32:02.9 

V37 

RRab 

17.947 

17.482 

0.83 

0.63 

0.519174 

- 

6399.4511 

0.5192 

16:46:58.87 

+47:32:11.9 

V38 

RRab-Bl 

18.06 

17.68 

1.06 

0.78 

0.522159 

- 

5665.3610 

0.5222 

16:46:57.09 

+47:31:11.4 

V39 

RRc 

18.113 

17.678 

0.42 

0.31 

0.332997 

- 

5725.3982 

0.4993 

16:47:03.05 

+47:31:08.1 

V40 

RRab-Bl 

17.58 

17.51 

0.51 

0.50 

0.591408 

- 

6399.2877 

0.5914 

16:47:02.37 

+47:31:58.8 

V41 

RRab-Bl? 

18.12 

17.51 

0.42 

0.31 

0.634608 

- 

5725.1939 

0.6347 

16:47:02.78 

+47:31:30.4 

V42 

RRab 

18.07 

17.54 

0.50 

0.37 

0.621730 

- 

5664.4199 

0.6218 

16:47:01.30 

+47:32:00.9 

V43 

RRab-Bl 

17.94 

17.41 

1.14 

0.73 

0.567649 

- 

5725.1641 

0.5677 

16:46:56.99 

+47:32:01.8 

V44 

RRc 

17.988 

17.598 

0.45 

0.33 

0.357365 

- 

6356.4737 

0.3573 

16:47:00.57 

+47:30:51.6 

V45 

RRab 

18.027 

17.486 

0.48 

0.33 

0.640061 

- 

6356.3490 

0.6401 

16:46:59.18 

+47:30:22.2 

V46 

RRab 

17.965 

17.433 

0.28 

0.22 

0.638741 

- 

6400.3493 

0.6388 

16:46:53.97 

+47:31:24.8 

V47 

RRc-Bl 

17.80 

18.03 

0.36 

0.40 

0.333948 

- 

5724.3047 

0.3339 

16:47:03.40 

+47:31:26.5 

V48 

RRc 

17.85 

17.47 

0.38 

0.26 

0.340270 

- 

5753.1385 

0.5164 

16:47:01.46 

+47:31:20.0 

V49 

RRab 

18.08 

17.51 

0.35 

0.24 

0.637320 

- 

5666.3954 

- 

16:46:54.75 

+47:32:36.5 

V50 

RRab-Bl 

18.11 

17.80 

1.25 

0.96 

0.492440 

- 

5752.3148 

- 

16:46:56.78 

+47:31:43.9 

V51 

RRab 

16.66 

15.65 

0.26 

0.11 

0.548110 

- 

5725.4078 

- 

16:46:57.21 

+47:31:27.4 

V52 

RRab 

18.13 

17.53 

1.05 

0.57 

0.491680 

- 

6356.3733 

- 

16:46:57.98 

+47:32:02.3 

V53 

RRab-Bl? 

18.12 

17.53 

0.36 

0.25 

0.625180 

- 

5753.2588 

- 

16:46:58.62 

+47:33:41.4 

V54 

RRab 

17.90 

17.34 

0.73 

0.47 

0.570370 

- 

6399.2877 

- 

16:46:59.20 

+47:31:50.2 

V55 

RRab 

18.00 

17.37 

1.19 

0.64 

0.525180 

- 

6462.3893 

- 

16:46:59.48 

+47:31:18.1 

V56 

RRab 

18.04 

17.34 

0.41 

0.28 

0.638310 

- 

6356.3490 

- 

16:46:59.70 

+47:31:52.6 

V57 

RRab 

18.22 

17.66 

0.99 

0.70 

0.561640 

- 

5753.1692 

- 

16:47:00.25 

+47:31:31.2 

V58 

RRab 

17.319 

16.527 

0.30 

0.16 

0.639500 

- 

5725.3496 

- 

16:47:00.34 

+47:32:14.6 

V59 

RRc 

18.0 

17.7 

0.55 

0.21 

0.335780 

- 

5663.4655 

- 

16:46:34.04 

+47:32:26.0 

V60 

RRc 

18.016 

17.684 

0.38 

0.25 

0.260380 

- 

6358.4103 

- 

16:46:55.27 

+47:31:31.0 

V61 

RRc-Bl 

17.67 

17.16 

0.31 

0.20 

0.343050 

- 

5753.2932 

- 

16:46:56.59 

+47:31:41.6 

V62 

RRc-Bl 

17.48 

16.81 

0.27 

0.16 

0.299040 

- 

6357.4281 

- 

16:46:59.52 

+47:32:02.0 

V63 

RRc 

16.750 

15.893 

0.12 

0.05 

0.341660 

- 

5752.2898 

- 

16:47:00.41 

+47:31:27.6 
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Table 3. Continued 


Variable 

Star ID 

Variable 

Type 

< V > 
(mag) 

< / > 
(mag) 

Av 

(mag) 

Aj 

(mag) 

P (days) 
this work 

p 

(d Myr-l) 

YUYimax 
(+2 450 000) 

P (days) 
(BCVOl) 

RA 

(J2000.0) 

Dec. 

(J2000.0) 

V64 

SR 

15.42 

13.94 

0.29 

0.13 

- 

- 

- 

- 

16:46:55.87 

-r47:31:59.9 

V65 

SR 

15.22 

13.65 

0.17 

0.12 

- 

- 

- 

- 

16:46:58.11 

+47:31:20.5 

V66 

SR 

15.04 

13.54 

0.33 

0.16 

- 

- 

- 

- 

16:46:59.01 

+47:31:58:0 

V67 

CWB 

17.12 

16.53 

1.35 

0.84 

1.51547 

- 

5724.1700 

- 

16:47:00.83 

+47:31:41.3 

V68 

SX Phe 

20.49 

20.48 

0.4 

- 

0.0384563 

- 

5663.3201 

- 

16:46:52.57 

+47:33:23.1 

V69 

7 

20.46 

20.24 

0.22 

- 

0.271264 

- 

6356.3416 

- 

16:47:01.85 

+47:32:13.3 

V70 

SR 

15.18 

13.55 

0.36 

0.21 

- 

- 

- 

- 

16:46:54.22 

+47:33:37.2 

V71 

SR 

15.22 

13.35 

0.39 

0.20 

- 

- 

- 

- 

16:46:58.58 

+47:31:36.9 

V72 

SR 

15.35 

13.99 

0.12 

0.10 

- 

- 

- 

- 

16:47:02.81 

+47:30:21.2 

V73 

CW? 

16.62 

15.78 

- 

- 

- 

- 

- 

- 

16:47:00.85 

+47:31:43.4 



X-pixel 

Figure 4. Minimum value of the string-length parameter Sq calculated for 
the 2113 stars with a light curve in our V reference image, versus the CCD 
X-coordinate. Colours are as in Fig.|^ The horizontal line is an arbitrarily 
defined threshold below which most variables seem to fall. All stars below 
the threshold were individually explored for variability. The only true vari¬ 
ables above the threshold are one SX Phe (V68), one unclassified variable 
blue straggler (V69), and one noisy RRab (V30). 

average heliocentric Julian day of the five images used to form the 
reference image, 2455665.92 d. 

3 VARIABLE STARS IN NGC 6229 
3.1 Search for new variable stars 

All previously known 48 variables were identified by their equato¬ 
rial coordinates listed in the CVSGC. The main approach to iden¬ 
tify new variables was the string-length method (Burke, Rolland & 
Boy 1970; Dworetsky 1983), in which each light curve is phased 
with periods between 0.02 and 1.7 day and a normalised string- 
length statistic Sq is calculated for each trial period. Extending the 
period base to values ~ 2.5d does not produce a significantly dif¬ 
ferent diagram. On the other hand if the period range is stretched to 
much longer periods, i.e. > lOd, the routine tends to find multiples 


of the true periods for some RR Lyrae. The long period variables 
are to be spotted by the other methods described below. In Fig. 
we plot the minimum Sq value for each light curve as a function 
of their corresponding CCD X-coordinate. The known variables 
are plotted with the coloured symbols described in the caption. The 
horizontal line is a threshold set by eye at Sq = 0.30, below which 
most variables seem to fall and then it is below this line where we 
might expect to find previously undetected variables. Then we care¬ 
fully inspected for variability the light curve of each star with Sq 
below 0.30. This approach allowed us to detect the new variables 
now labelled V49-V67 and to confirm the non-variability of stars 
V26 and V28. Despite falling below the threshold, V25 does not 
show convincing signs of variability. We also note that this method 
did not work for the SX Phe star, a fact already pointed out before 
by Arellano Ferro et al. (2010) in the study of NGC 5024. 

A second approach towards searching for variable stars was 
to separate the light curves of stars in a given region of the CMD 
where variables are expected, e.g. HB, blue stragglers region, upper 
instability strip and the RGB. A further more detailed inspection of 
stars in these regions revealed two new variables in the blue strag¬ 
glers region, V68 and V69. At least star V68 can be classified as 
SX Phe. Also the SR variables labelled V70-V72 and a possible 
CW star V73, which we had missed from the string-length method, 
were found in this way. 

Finally, we used a third approach that consisted in detecting 
PSF-like peaks in a stacked image built from the sum of the abso¬ 
lute valued difference images normalised by the standard deviation 
in each pixel as described by Bramich et al. (2011). This method 
allowed us to confirm the variability of all of the new variables dis¬ 
covered by the previous methods but no new variables emerged. 

The confirmed variables, their mean magnitudes, amplitudes, 
periods, times of maximum light and equatorial coordinates are 
listed in Table Detailed finding charts of all variables are given 
in Figs. and The light curves of RRL stars, phased with the 
ephemerides given in Tableare displayed in Figs|7]and[^ Com¬ 
ments on individual stars are made in Appendix A. 

From the above methods of variable star finding we estimate 
that our search is complete down to about V ~ 18.5 for amplitudes 
larger than 0.03 mag and periods between 0.2 d and a few tens 
of days in the outer regions of the cluster. However, we note that 
the core is very concentrated and our subtractions there are noisy. 
Therefore, our variable search in the core is not complete. In two 
recent papers Skottfelt et al. (2013; 2015) have shown the benefits 
of using electron multiplying CCD (EMCCD) photometry and DIA 
to extract faint variables in the dense central regions of globular 
clusters. This is achieved from a high frame-rate time-series and by 
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Figure 5. Finding charts constructed from our V reference image; north is up and east is to the right. The cluster images are 7.3 x4.4 and 1.25x0.90 arcmin^. 


shifting and adding frames to counteract the blurring atmospheric 
effects which contributes towards a high spatial resolution. Thus 
we believe the core of NGC 6229 would be an excellent target for 
EMCCD photometry. 


3.2 RR Lyrae stars with a secular period variation 

The combination of our data with those from Baade (1945) and 
Mannino (1960), available for variables V1-V21, give a time-base 
of about 81 years. For the benefit of future work on these variables 
we have uploaded these archival data into the Strasbourg astronom¬ 
ical Data Center (CDS|^ For these stars the light curves of the 
archival data phased with the ephemerides in Tablej^are also shown 

^ data are available at the following links: http://cdsarc.u-strasbg.fr/cgi- 
bin/VizieR?- source=1/ Ap J/102/17 and http: //cdsarc. u- strasbg. fr/cgi- 
bin/VizieR?-source=J/other/MmS AI/31.187 


in Figs.[7]and[^ Prominent phase shifts are visible for the different 
data sets in several variables, suggesting secular period variations. 
Despite the large time-base, the number of times of maximum that 
can accurately be determined for each variable is limited, prevent¬ 
ing us from calculating the secular period change with the classical 
0-C residuals analysis. Alternatively, we used the approach intro¬ 
duced by Kains et al. (2015) for the case of NGC 4590, which for 
completeness is described below. The period change can be repre¬ 
sented as 

P(t) = Po + ^(t-^), (5) 

where /3 is the period change rate expressed in units of d d“^, and 
Po is the period at the epoch E. The number of cycles at time t 
elapsed since the epoch E is: 
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hence the phase at time t can be calculated as 

(l,{t) = NE-[NE\. (7) 

Then we construct a grid (P,/3) with values in steps of 10(7 
and 0.004 d Myr“^ in P and [3 respectively, around a given starting 
period Po, and examine the light curve dispersions for each pair. We 
choose as a solution the one that produces the minimum dispersion 
of our data combined with those of Baade (1945) and Mannino 
(1960). In Fig. we show three examples of the light curves phased 
with eq.[^ those of the RRab stars VI and V2 and the RRc V15. 
The resulting P and /3 for those stars with an evident period change 
are listed in Tablej^ In a few cases where the archival data are very 
scattered, we have not been able to find a reliable solution (e.g. for 
stars VI7 and VI8); however, their secular variation is clear and 
the proper calculation will have to wait until sufficient data become 
available. 

The period change rates that we found are very scattered; they 
average —0.116di0.201d Myr“^. While this number is consistent 


with zero, the scatter is very large. Taking the sign of the period 
change rate as an indication of the direction in which the star is 
evolving, we could not discern for NGC 6229 a systematic direc¬ 
tion towards which stars are evolving. A summary of the observed 
period changes in RR Lyrae stars in globular clusters is offered by 
Catelan (2009), in particular his Fig. 15, according to which the pe¬ 
riod change rate (5 is related to the £ = (B — R)/(B + V + R) 
Lee-Zinn parameter, defined to describe the stellar distribution on 
the HB (Zinn 1986; Lee 1990), where B, V, R refer to the number 
of stars to the blue, inside and to the red of the RR Lyrae instability 
strip. Catelan’s diagram shows that, for clusters with C < 0.5 the 
value of /3 is negligible, and it increases exponentially for larger £, 
i.e. for clusters with blue MB’s. From star counts in our CMD we 
estimate for NGC 6229 C — 0.18, which can be compared with the 
value 0.24 reported by Mackey & van den Bergh (2005), and that 
would imply ^ ~ 0. Despite this the average /S found in this paper 
is very negative but with a large scatter. This may imply that rather 
erratic period changes dominate the period changes produced oth¬ 
erwise by stellar evolution across the instability strip in NGC 6229. 
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Figure 7. Light curves of the RRab stars in NGC 6229. They are phased with the ephemerides listed in Table the periods in days are given next to the 
identification. To preserve the vertical scale and limits the vertical axes display V— < V > and I— < I >; the individual mean values are also listed in the 
table. Orange circles represent the data from Baade (1945) and green circles are data from Mannino (1960), whenever available. Phase differences between 
the archival data and our data suggest secular period variations that are calculated in § |3.2| Variables with names V49 and onwards are new discoveries in this 
paper. 
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Figure 7. Continued 
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Figure 8. Light curves of the RRc stars in NGC 6229 phased with the periods listed in Table[^ 
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Figure 9. Examples of light curves phased with the ephemerides including the period change rate (3 which 
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3.3 Bailey diagram and Oosterhoff type 

The periods of the 42 RRab and 15 RRc stars listed in Tablej^aver- 
age 0.561d and 0.322d, respectively. These values clearly indicate 
that NGC 6229 is an Oosterhoff type I cluster. 

The Bailey diagram for the V and I filters is shown in Fig.p^ 
The distribution of the RRab stars in each plane is in close agree¬ 
ment with the loci identified for Oo I clusters by Cacciari et al. 
(2005) (M3), Kunder et al. (2013a) (NGC 2808) and Arellano Ferro 
et al. (2014) (NGC 3201). A number of outliers are noted: the case 
of V51 can be understood since its amplitude is uncertain due to its 
Blazhko nature. For V14 and V37, with clean non-modulated light 
curves, we do not have an explanation. In the upper end we note 
VI is the only non-modulated star sitting near the locus for evolved 
stars, however its position on the HB (see the bottom panel of Fig. 

does not add to that interpretation. The Bailey diagram further 
confirms NGC 6229 as an Oo I type cluster. 

The distribution of RRc stars is expected to have a parabolic 
form (Bono et al. 1997) as in fact it has been found in some Oo 
II clusters e.g. NGC 2419 (Di Criscienzo et al. 2011), M9 (Arel¬ 
lano Ferro et al. 2013), M55 (Olech et al. 1999), M30 (Kains et al. 
2013) and Kunder et al. (2013b) for a group of 14 Oo II clusters. 
We have assembled a sample RRc stars in Oo I clusters (M3, M5, 
M79, NGC 6229 and NGC 6934) with determined amplitudes in V 
and/or I and without obvious amplitude modulatios. From their dis¬ 
tribution (not shown here) we calculated the red parabolas shown 
in Fig.[^ The scatter about these parabolas is sometimes consid¬ 
erable, like in the case of Oo II clusters of Kunder et al. (2013a). 
The mathematical expresions for these parabolas are: 

=-3.95 + 30.17P-51.35P^ (8) 


=-2.72 + 20.78P-35.51P7 (9) 


4 RR LYRAE STARS: [FE/H] AND Mv FROM LIGHT 
CURVE FOURIER DECOMPOSITION 

Stellar physical parameters, such as [Fe/H], My, Tgff, mass and 
radius for RR Lyrae stars can be calculated via the Fourier decom¬ 
position of their V light curves into its harmonics as: 


m{t) = Ao + cos (^ k {t- E) + 0/e), (10) 

k=i 

where m{t) are magnitudes at time t, P is the period and E the 
epoch. A linear minimization routine is used to derive the ampli¬ 
tudes Ak and phases (pk of each harmonic, from which the Fourier 
parameters cpij = j(pi — i<pj and Rij — AijAj are calculated. The 
mean magnitudes Aq, and the Fourier light curve fitting parameters 
of individual RRab and RRc stars are listed in Table |4] In this ta¬ 
ble we have excluded stars with evident amplitude modulations or 
excessive noise. 

These Fourier parameters and the semi-empirical calibra¬ 
tions of Jurcsik & Kovacs (1996), for RRab stars, and Morgan 
et al. (2007), for RRc stars, are used to obtain [Fe/H]zw on the 
Zinn & West (1984) metallicity scale which have been trans¬ 
formed to the UVES scale using the equation [Fe/H]uvES= —0.413 
+0.130[Fe/H]zw - 0.356[Fe/H]|w (Carretta et al. 2009). The 
absolute magnitude My can be derived from the calibrations of 
Kovacs & Walker (2001) for RRab stars and of Kovacs (1998) for 
the RRc stars; these two calibrations are naturally based on dif¬ 
ferent sets of calibrators in assorted Galactic and LMC globular 
clusters and are therefore independent. However, their zero points 
have been tied to the LMC distance modulus of 18.5zb0.1 mag (e.g. 
Clementini et al. 2003) as discussed in detail by Arellano Ferro 
et al. (2010). The effective temperature Tgff was estimated using 
the calibration of Jurcsik (1998). For brevity we do not explicitly 
present here the above mentioned calibrations; however, the corre¬ 
sponding equations, and most importantly their zero points, have 
been discussed in detail in previous papers (e.g. Arellano Ferro et 
al. 2011; 2013) and the interested reader is referred to them. 

It is pertinent to recall that the calibration for [Fe/H] for RRab 
stars of Jurcsik & Kovacs (1996) is applicable to RRab stars with 
a deviation parameter Dm, defined by Jurcsik & Kovacs (1996) 
and Kovacs & Kanbur (1998), not exceeding an upper limit. These 
authors suggest Dm < 3.0. The Dm is listed in column 10 of 
Table |4] However we have relaxed the criterion to include three 
stars with 3.7 < Dm <4.1. Those RR Lyrae that fall off the HB, 
discussed in §|^ are not included for the luminosity calculation as 
indicated in the footnote of Table [5] 

It is noteworthy at this stage that there is a new calibration of 
[Fe/H]UVES for RRc stars, calculated for an extended number of 
stars and clusters by Morgan (2013). As described in her paper the 
results of the new calibration are in general consistent with the old 
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Table 4. Fourier coefficients for = 0,1, 2, 3, 4, and phases 02i, 03i and 041, for RRab and RRc stars. The numbers in parentheses indicate the 
uncertainty on the last decimal place. Also listed is the deviation parameter (see Section]^. 


Variable ID Aq 

(y mag) 

Al 

(V mag) 

A 2 

(V mag) 

A 3 

(V mag) 

A 4 

(V mag) 

021 

031 

041 


RRab stars 

VI 

17.936(1) 

0.383(2) 

0.206(2) 

0.132(2) 

0.088(2) 

3.988(12) 

8.316(18) 

6.422(26) 

0.9 

V5 

18.160(2) 

0.359(2) 

0.174(2) 

0.128(2) 

0.081(2) 

3.876(18) 

8.099(27) 

6.153(37) 

2.0 

VIO 

18.101(2) 

0.332(2) 

0.158(2) 

0.120(2) 

0.075(2) 

3.894(18) 

8.198(25) 

6.232(38) 

1.0 

Vll 

18.119(3) 

0.302(3) 

0.151(3) 

0.102(2) 

0.065(3) 

4.043(29) 

8.386(43) 

6.526(63) 

0.4 

V13 

18.106(2) 

0.308(3) 

0.159(3) 

0.112(2) 

0.069(2) 

3.886(25) 

8.238(36) 

6.235(53) 

1.4 

V14 

17.840(2) 

0.366(3) 

0.182(3) 

0.125(3) 

0.078(3) 

3.737(22) 

7.790(32) 

5.681(48) 

1.9 

V18 

18.066(2) 

0.273(2) 

0.128(2) 

0.091(2) 

0.052(2) 

4.124(24) 

8.420(36) 

6.823(56) 

3.7 

V19 

18.130(2) 

0.432(3) 

0.209(3) 

0.153(3) 

0.101(3) 

3.811(20) 

7.957(29) 

5.906(42) 

1.5 

V21 

18.063(2) 

0.304(3) 

0.152(3) 

0.105(3) 

0.070(3) 

4.011(32) 

8.293(47) 

6.308(95) 

1.5 

V35 

18.062(1) 

0.115(2) 

0.039(2) 

0.020(2) 

0.002(2) 

4.259(64) 

8.924(111) 

7.049(84) 

6.0 

V37 

17.947(2) 

0.282(3) 

0.138(3) 

0.095(3) 

0.057(3) 

3.959(29) 

8.171(43) 

6.237(63) 

1.1 

V45 

18.027(1) 

0.174(2) 

0.073(2) 

0.037(2) 

0.014(2) 

4.190(31) 

8.794(55) 

7.540(13) 

26.9 

V46 

17.965(1) 

0.115(2) 

0.039(2) 

0.018(2) 

0.009(2) 

4.280(31) 

9.042(55) 

7.899(13) 

32.5 

V51 

16.658(1) 

0.068(2) 

0.040(2) 

0.031(2) 

0.022(2) 

4.050(71) 

8.312(101) 

6.436(183) 

4.1 

V58 

17.319(1) 

0.102(1) 

0.046(1) 

0.023(1) 

0.007(1) 

4.046(14) 

8.664(14) 

7.056(13) 

4.1 

RRc stars 

V15 

18.118(2) 

0.252(2) 

0.056(2) 

0.015(2) 

0.014(2) 

4.635(43) 

2.716(158) 

1.449(167) 

- 

V23 

17.868(1) 

0.185(2) 

0.009(2) 

0.011(2) 

0.007(2) 

5.742(211) 

4.718(183) 

3.143(270) 

- 

V24 

18.016(1) 

0.217(2) 

0.030(2) 

0.015(2) 

0.014(2) 

4.522(64) 

3.262(124) 

2.189(134) 

- 

V36 

17.942(2) 

0.149(2) 

0.005(2) 

0.002(2) 

0.003(2) 

4.609(464) 

4.956(1296) 

3.515(755) 

- 

V39 

18.113(2) 

0.210(2) 

0.032(2) 

0.018(2) 

0.011(2) 

4.695(71) 

3.834(127) 

2.525(209) 

- 

V44 

17.988(1) 

0.216(2) 

0.027(2) 

0.008(2) 

0.010(2) 

4.883(73) 

4.051(246) 

2.750(185) 

- 

V48 

17.849(2) 

0.182(3) 

0.019(2) 

0.009(2) 

0.006(2) 

5.726(130) 

4.425(269) 

1.978(420) 

- 

V60 

18.016(2) 

0.176(2) 

0.027(3) 

0.005(2) 

0.003(3) 

4.402(95) 

2.649(481) 

0.893(751) 

- 

V63 

16.750(1) 

0.046(1) 

0.002(1) 

0.004(1) 

0.004(1) 

4.550(89) 

4.105(328) 

2.100(362) 

- 


calibration from Morgan et al. (2007), and that the average differ¬ 
ences are 0.02 dex. We have also carried out the calculation using 
the new calibration and found [Fe/H]uvES = —1.37, i.e. 0.08 dex 
more deficient than the value from the calibration of Morgan et al. 
(2007), which is in better agreement with the result from the RRab 
calibration. 

The physical parameters for the RR Lyrae stars and the 
inverse-variance weighted means are reported in Tablej^ Two inde¬ 
pendent estimations of [Fe/H]uvES have been found from 12 RRab 
and 8 RRc stars in NGC 6229: -1.31 ± 0.01 and -1.29 ± 0.12 
respectively. 

The weighted mean My values for the RRab and RRc stars 
are 0.621 ±0.001 mag and 0.581 ±0.001 mag respectively and will 
be used in section to estimate the mean distance to the parent 
cluster. 


5 ON THE STRUCTURE OE THE HORIZONTAL 
BRANCH 

In this section we analyse the distribution of RRL stars in the HB. 
The expanded HB in Fig. [^displays in detail the positions of most 
of the RR Lyrae stars. RRab and RRc are distinguished as well as 
Blazhko and stable variables. There is a neat segregation between 
stable RRab and RRc with a border at (U — /)o ~ 0.45. Similarly 
clear segregations have been found before in NGC 5024 (Arellano 
Ferro et al. 2011; their Fig. 4) and NGC 4590 (Kains et al. 2015; 
their Fig. 11) both at (U-/)o ~ 0.45—0.46. In the bottom panel of 
Fig.j^the corresponding border lines for NGC 5024 and NGC 4590 


duly reddened are indicated. It can be seen that the border line be¬ 
tween stable RRab and RRc is the same in the three clusters. We 
should emphasise that the above empirical border is trespassed by 
a few Blazhko variables (NGC 6229) and the double mode RRL 
(NGC 4590). Looking for this segregation effect in clusters previ¬ 
ously studied by our group, we checked the CMDs and noted that 
with very few exceptions the effect is also visible in NGC 5053 
(Arellano Ferro et al. 2010; Fig. 8) (except for VIO with a poor I 
light curve), NGC 1904 (Kains et al. 2012) and NGC 7099 (Kains 
et al. 2013), despite their low number of RRL. On the other hand 
the RRab-RRc splitting is definitely not seen in NGC 3201 (Arel¬ 
lano Ferro et al. 2014). 

It is illustrative to display the HB in the Teff-logL plane. In or¬ 
der to convert the observational HB shown in Fig.j^into the HR di¬ 
agram, we have proceeded as in Arellano Ferro et al. (2010) where 
a polynomial for Teff-(U — /)o is given based on the HB models of 
VandenBerg et al. (2006) with the colour-Tgff relations as described 
by VandenBerg & Clem (2003). Fig. illustrates the distribution 
of RRL using the same colour code as in Fig.[^ The red and blue 
borders of the instability strip (IS) for the fundamental mode (F) 
and first overtone (FO), calculated by Nemec et al. (2011) based on 
the Warsaw pulsation code for a mass of 0.65 M©, are shown in 
the figure. The two vertical black lines correspond to the empirical 
border between RRab and RRc stars seen in NGC 5024 (solid) and 
NGC 4590 (dashed) which match the evident border in NGC 6229 
and that, in our opinion, represent the empirical red border of the 
first overtone instability strip (FORE). This is some 100-200K hot¬ 
ter than the theoretical border found by Nemec et al. (2011). An- 
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Table 5. Physical parameters for the RRab and RRc stars. The numbers in parentheses indicate the uncertainty on the last decimal place. 


Star 

[Fe/H]uvES 

Mv 

log T'eff 

log(L/L 0 ) 

M/M© 

R/Rq 

RRab stars 

VI 

-1.385(18) 

0.511(3) 

3.809(7) 

1.696(1) 

0.70(6) 

5.69(1) 

V5 

-1.394(27) 

0.611(3) 

3.812(8) 

1.656(1) 

0.70(7) 

5.36(1) 

VIO 

-1.378(25) 

0.604(3) 

3.810(8) 

1.658(1) 

0.68(7) 

5.42(1) 

Vll 

-1.215(39) 

0.612(4) 

3.811(10) 

1.655(2) 

0 . 66 ( 8 ) 

5.40(2) 

V13 

-1.302(34) 

0.636(4) 

3.811(9) 

1.645(2) 

0.66(7) 

5.32(1) 

V14 

-1.437(32) 

0.711(4) 

3.817(9) 

1.616(2) 

0.71(7) 

5.00(1) 

V18 

-1.245^3) 

0.614(3) 

3.806(9) 

1.655(1) 

0.66(7) 

5.50(1) 

V19 

-1.297(27) 

0.640(4) 

3.820(8) 

1.644(2) 

0.73(7) 

5.10(1) 

V21 

-1.317^5) 

0.610(4) 

3.810(13) 

1.656(2) 

0 . 66 ( 10 ) 

5.42(2) 

V35“ 

- 

0.664(3) 

3.804(12) 

1.634(1) 

0.56(8) 

5.43(1) 

V37 

-1.254(40) 

0.696(4) 

3.813(10) 

1.622(2) 

0.65(8) 

5.15(2) 

V45“ 

- 

0.602(3) 

3.799(7) 

1.659(1) 

0.63(5) 

5.72(1) 

V46“ 

- 

0.657(3) 

3.800(21) 

1.637(1) 

0.58(15) 

5.55(1) 

V51^ 

-1.229(92) 

- 

3.804(21) 

- 

- 

- 

V58^ 

-1.243(13) 

- 

3.797(7 ) 

- 

- 

- 

Weighted 

Mean 

-1.31(1) 

0.621(1) 

3.810(2) 

1.652(1) 

0 . 66 ( 2 ) 

5.43(1) 


RRc stars 


V15 

-1.22(25) 

0.603(2) 

3.872(1) 

1.659(1) 

0.61(1) 

4.09(1) 

V23 

-1.49(47) 

0.464(9) 

3.859(1) 

1.714(4) 

0.47(1) 

4.62(2) 

V24 


0.577(3) 

3.868(1) 

1.669(1) 

0.56(1) 

4.21(1) 

V36“ 

- 

0.660(21) 

3.887(8) 

1.636(8) 

0.50(5) 

3.72(4) 

V39 

-1.34(26) 

0.555(3) 

3.865(1) 

1.678(1) 

0.52(1) 

4.31(1) 

V44 

-1.50(58) 

0.527(3) 

3.862(1) 

1.689(1) 

0.50(1) 

4.43(1) 

V48 

-1.03 be) 

0.524(6) 

3.867(2) 

1.690(2) 

0.51(1) 

4.33(1) 

V60 

-1.11(69) 

0.673(4) 

3.874(3) 

1.631(2) 

0.58(2) 

3.92(1) 

V63^ 

-1.27^6) 

- 

3.865(2) 

- 

- 

- 

Weighted 

Mean 

-1.29(12) 

0.581(1) 

3.867(1) 

1 . 668 ( 1 ) 

0.54(1) 

4.18(1) 


No [Fe/H] estimate since Dm > 4.1. V36 has a very peculiar value of (psi; No My estimate since star falls off the HB. 


other feature is that the fundamental red edge (FRE) in NGC 6229 
extends further to the red by at least 400K. 

Regarding the distribution of RRc stars, they are well con¬ 
strained to the first overtone instability strip. The distribution of 
RRL on the HB observed in NGC 6229 can be explained by the ar¬ 
guments of Caputo et al. (1978) which involve the occurrence of a 
hysteresis mechanism (van Albada & Baker 1973) for stars crossing 
the IS. According to this mechanism the stars in the ”either-or” re¬ 
gion can retain the mode they were pulsating in before entering the 
region which depends on whether the star is coming from the blue 
side as a RRc or from the red side as a RRab. Furthermore, Caputo 
et al. (1978) suggest that the original mode in which a RRL is pul¬ 
sating depends on the location of the starting point at the Zero Age 
Horizontal Branch (ZAHB), which in turn depends on the mass and 
the chemical composition of the star; they used these ideas to ex¬ 
plain the existence of the two Oosterhoff groups. A schematic and 
clear explanation of this can be found in their Fig. 3. In summary 
there are two scenarios; 1) the ZAHB point is in the fundamental 
zone, leading to an assortment of RRc and RRab in the ”either-or” 
region and a lower value of the average period for both RRc and 
RRab and a lower proportion of RRc stars, hence an Ool cluster, 
and 2) the ZAHB point is bluer than the fundamental zone, in the 
”either-or” or in the first overtone regions, then the ”either-or” re¬ 
gion is populated exclusively by RRc stars and the RRab stars are 


to be found only in the fundamental region, producing larger av¬ 
erages of periods and a higher proportion of RRc and hence OoII 
type cluster. 

The above scenario seems to explain well the clean segrega¬ 
tion of RRc and RRab in OoII clusters with evolution across the 
IS towards the red. In agreement with this picture we can mention 
the analysis of Pritzl et al. (2002b). These authors have shown how 
in clusters with blue HBs, i.e. large values of £, stars with masses 
below a critical mass on the ZAHB to the blue of the IS, evolve 
redwards and spend sufficient time to contribute to the population 
of RRL. These arguments suggest then that clusters with a blue 
HB morphology are OoII, with redwards evolution, which in turn 
favours the segregation between RRc and RRab at the IS. 

We have explored the CMD of several OoII clusters published 
by our group over the last few years (e.g. NGC 288, NGC 1904, 
NGC 5024, NGC 5053, NGC 5466, NGC 6333, NGC 7099) and 
noted the clear segregations, with a few exceptions of indiviual 
stars with noisy light curves. The one exception is NGC 4590, a 
clear OoII with a red HB (£=0.17) and evidences, from the secular 
period changes analysis of 14 RRL, that evolution does not occur 
in a preferential direction (Kains et al. 2015). On the other hand 
the Ool cluster NGC 3201 (£=0.08) as expected, shows no signs 
of segregation (Arellano Ferro et al. 2014). Likewise NGC 6229 
has a low £ (0.24) and no preferential evolution direction from the 
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log P 

Figure 10. The log P vs. amplitude plane for the V and / amplitudes of the 
RRL in NGC 6229. Filled and open symbols represent RRab and RRc stars 
respectively. Triangles are used for stars with Blazhko modulations. In the 
top panel the continuous and segmented lines are the loci found by Cacciari 
et al. (2005) for unevolved and evolved RRab stars, respectively, in the Ool 
cluster M3. The black parabola was found by Kunder et al. (2013b) from 14 
Oo II clusters. We calculated the red parabola from a sample of RRc stars 
in hve Oo I clusters and avoiding Blazhko variables. In the bottom panel 
the black segmented locus was found by Arellano Ferro et al. (2011, 2013) 
for the OoII clusters NGC 5024 and NGC 6333. The blue loci are from 
Kunder et al. (2013a) for the Ool cluster NGC 2808. The red parabola was 
calculated from RRc stars in three Oo I clusters with / amplitudes. 

scattered period change rate found in § |3.2| in spite of this the RRc 
and the RRab are neatly separated. 

A close inspection of the distribution of clusters of all Ooster- 
hoff types on the [Fe/H]-£ plane from the data in Table 3 of Cate- 
lan (2009), reveals NGC 4590 as the OoII cluster with the reddest 
HB with £ < 0.4 and that there are five Ool clusters (NGC 6266, 
NGC 6284, NGC 6402, NGC 6558 and NGC 6626) with blue HB’s 
(£ > 0.4). It would be of great interest to explore their RRL distri¬ 
butions in the IS. 

One more point of interest in Fig. m is the distribution of 
Blazhko variables. We have plotted with triangles only stars that 
show clear Blazhko amplitude modulations while marginal cases 
(e.g. V50, V52, V54) are plotted as stable stars. In his paper Gillet 
(2013) concluded that the Blazhko effect is not a consequence of 
the transition between modes and that a star that penetrates the 



Figure 11. The HR diagram of the horizontal branch of NGC 6229. The 
coloured lines are the IS borders for the fundamental mode (FM) and the 
first overtone (FO) according to Nemec et al. (2011). The two vertical black 
lines correspond to the empirical border between RRab and RRc stars seen 
in NGC 5024 (solid) and NGC 4590 (dashed). See also Fig.j^and the text 
for a discussion. 

’’either or” region (RRc or RRab) becomes a Blazhko variable, as 
seems to be supported by the field Blazhko RRL plotted in Gillet’s 
Fig. 1. The later statement implies that all stars in the ’’either 
or” region should be Blazhko variables, which is not observation- 
ally supported by the distribution of Blazhko stars in NGC 4590, 
NGC 5024 and NGC 6229. Some outlier Blazhko variables from 
the IS in NGC 6229, mostly RRab stars, are likely not the result 
of hysteresis in the change of mode process as once suggested by 
Arellano Ferro et al. (2010) but rather a consequence of multi¬ 
shock perturbation at the main pulsation instability yet to be ex¬ 
plained by theory (e.g. Gillet 2013). 

In the CMD of Fig. we have overplotted two versions of 
the ZAHB for [Fe/H]=-1.42 and [a/Fc]=+0A; for Y=0.25 and 
Y=0.27, calculated from the Victoria-Regina stellar models of Van- 
denBerg et al. (2014). The model for Y=0.25 is a reasonable lower 
envelope to the distribution of RRL. The difference with a non- 
diffusive model for [Fe/H]=—1.31, [a/Fe]=-i-0.3 and Y=0.25 from 
VandenBerg et al. (2006) is negligible. Apparently the change in 
the primordial helium content largely compensates for the differ¬ 
ence between diffusive and non-diffusive models (VandenBerg pri¬ 
vate communication). The only two stable RRab’s a little below this 
ZAHB are V33, which shows some scatter near minimum, and V55 
whose flux is likely contaminated by its neighbour V31 (see the dis¬ 
cussion in § |A0.1| ). If the intrinsic thickness of the stellar distribu¬ 
tion on the HB is between 0.1-0.2 mag, e.g. the case of M3 in Fig. 
1 of Cacciari et al. (2005), it is clear that a number of RRL in NGC 
6229 are overluminous. A natural explanation would be evolution 
off the ZAHB, however we did not find evidence of stars evolved in 
the amplitude-period diagram of Fig. Also no evidence of sys¬ 
tematic secular period increases in this rather red-HB cluster have 
been found (see § |3.2| ). Furthermore, stars in an advanced evolu¬ 
tion towards the AGB should be brighter and have longer periods, 
however we find no correlation between <V> and period for the 
RRL in NGC 6229. 

There may be other reasons for overluminosity than evolution 
off the ZAHB. For example, helium abundance dispersion in HB 
stars could be produced if more than one generation of stars are 
mixed in a given GC. D’Antona & Caloi (2008) have argued that 
second-generation stars, formed from enriched material of first- 
generation progenitors, may have a spread in Y above the primor- 
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Figure 12. Period-colour diagram of stable RRL stars. Filled circles are 
used for RRab stars and empty circles for RRc stars. The periods of the 
RRc stars were fundamentalized by adding 0.127 to logP. The reduced 
period P' is defined as logP' = logP + 0.336(< V > —18.021). See 
text for further discussion. 


.N' 


• • 


V14 


dial value, and that these stars will occupy a bluer position on the 
HB than the first-generation ones, which sit to the red of the in¬ 
stability strip. It has been noted by Catelan et al. (2009) that HB 
stars with higher Y should also be brigther due to their stronger 
H-burning shells, and they have searched for such an effect in M3. 
By comparing with HB models with Y in the range 0.23-0.28, they 
concluded that no systematic overluminosity in blue HB stars is ev¬ 
ident, and hence the Y spread in M3 should not be larger than 0.01. 
For NGC 6229 we have overplotted ZAHB models from Vanden- 
Berg (2014) for Y = 0.25, 0.27, and 0.285, and found no evidence 
of the blue HB stars being brighter than the red ones, so no spread 
in the value of Y is immediately obvious in this cluster either. Sim¬ 
ilar tight constraints would apply to the case of helium mixing on 
the RGB phase (Sweigart 1997). 

Another reason may be the presence of an unseen compan¬ 
ion. In this context it is of interest to explore the period-colour 
relation for stable RRL of Fig. In this diagram the period 
logP' = logP + 0.336(< V > —18.021) where 18.021 mag 
is the average <V > for all the stable RRL. The concept of logP' 
was introduced by van Albada & Baker (1971) and it was later 
slightly modified by Bingham et al. (1984). It significantly reduces 
the scatter in the period-colour diagram and is useful for identify¬ 
ing stars with peculiar photometric properties such as colours too 
red or too blue, and luminosities too bright or too faint relative to 
what is expected from their periods. In Fig. we note some out¬ 
liers and have labelled them. V33 has been commented on above, 
while V36 was believed to be an RRd in the CVSGC but classi¬ 
fied here as an RRc with a very short period among RRc stars. It 
is worth noting that V14, V31, V54 and V55 are also the brightest 
stable RRab in the HB (see the bottom panel of Fig.[^. Their po¬ 
sition on Fig, implies that they have too short a period for their 
colour, which implies that they have larger gravity and should be of 
lower luminosity. The fact that instead they are among the brightest 
in the HB, implies that they have an extra source of luminosity and 
this may be an unseen companion. It was through the same lines 
of argumentation that Cacciari et al. (2005) seem to have identified 
three RRL in M3 (V48, V58 and V146) that may have companions. 

Thus it seems likely that not only evolution has contributed 
to the large dispersion of RRL in the HB of NGC 6229 but also a 
dispersion in the primordial helium abundance and in some cases 
the presence of unseen companions. 


6 NGC 6229 DISTANCE AND METALLICITY EROM ITS 
VARIABLE STARS 

Four independent estimates of the distance to the cluster are pos¬ 
sible from our data; from the weighted mean My, calculated for 
the RRab and RRc from the Fourier light curve decomposition (Ta¬ 
ble 1^, which can be considered as independent since they come 
from different empirical calibrations; from the one known SX Phe 
star and the P-L relation; and from the tip of the RGB bolometric 
magnitude calibration. Below we expand on these solutions. 

Since NGC 6229 is a high Galactic latitude object and de¬ 
spite its large distance, the interstellar extinction and reddening are 
negligible, we adopted E{B — I/)=0.01 (Harris 1996). Given the 
mean My for RRL in Table|^we found a true distance modulus of 
17.385 zb 0.117 mag and 17.433 zb 0.096 mag using the RRab and 
RRc stars respectively, which correspond to the distances 30.0zb 1.5 
and 30.7 zb 1.1 kpc. The quoted uncertainties are the standard devi¬ 
ation of the mean. The distance to NGC 6229 listed in the catalogue 
of Harris (1996) (2010 edition) is 30.5 kpc in good agreement with 
our calculations. 

A third independent estimate of the distance to the cluster is 
via the one SX Phe identified in this work, V68, assuming that its 
period 0.0384563d corresponds to the fundamental mode. We have 
used the P-L relation of Cohen & Sarajedini (2012); 


My = -(1.640 zb 0.110) - (3.389 zb 0.090)logP/ (11) 

where Pf refers to the fundamental mode given in days. The above 
equation implies a true distance modulus of 17.303zb0.168 or a dis¬ 
tance of 28.9 zb 2.2 kpc. The uncertainty was calculated from cal¬ 
ibration errors neglecting the uncertainty of the period. It is worth 
noting that eq. [m when applied to a sample of S Scuti stars in the 
LMC produces a distance modulus of 18.49zb0.10 (Cohen & Sara¬ 
jedini 2012) which is consistent with the modulus 18.5zb0.1 used to 
set the zero points of the absolute magnitude calibrations for RRab 
and RRc stars (see Arellano Ferro et al. 2010 § 4.2). Thus the dis¬ 
tances derived here for the RRL and the SX Phe are on a consistent 
distance scale. 

Yet another approach to the cluster distance determination is 
by using the tip of the RGB. This method was developed for es¬ 
timating distances to nearby galaxies (Lee, Freedman & Madore, 
1993). Here we used the calibration of Salaris & Cassisi (1997) 
for the bolometric magnitude of the tip of the RGB in terms of the 
overall metallicity: 


= -3.949 - 0.178[M/H] + 0.008[M/H]^ (12) 

where [M/H] = [Fe/H] + log(0.638f + 0.362) and log f = [a/Fe] 
(Salaris et al. 1993). 

We found that the method is extremely sensitive to the star 
selection as the tip of the RGB, hence we restricted our calculation 
to the five reddest stars in the red box in the CMD of Fig.[^ We have 
assumed [Q(/Fe]=-i-0.4 (e.g Jimenez & Padoan 1998) and found an 
average distance of 34.9 zb 2.4 kpc. 

Within their own uncertainties the above distance determina¬ 
tions are satisfactorily in agreement. However, since the SX Phe 
determination is based on only one star and the tip of the RGB 
approach is very sensitive to the selected RGB tip stars, which in 
turn is scattered, we believe that the distance determination from 
the RRL stars is the most accurate. 

The average value of [Fe/H]uvES = — 1.31zb0.12 found from 
the Fourier decomposition of the light curves of 12 RRab and 8 RRc 
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stars in NGC 6229 compares very well with the numerous determi¬ 
nations in the literature. In the present edition of his catalogue of 
globular clusters Harris (1996) lists an average [Fe/H]zw = —1.47 
(or [Fe/H]uvES = -1.37) for NGC 6229. A list of individual re¬ 
sults from an assortment of indicators can be found in Table 7 of 
Borissova et al. (1999); they range from —1.26 to —1.44 in the ZW 
scale or from — 1.14 to — 1.34 in the UVES scale with uncertainties 
running within 0.08-0.15 dex. 


6.1 A brief comment on the age of NGC 6229 

We have not attempted an independent estimation of the age of 
NGC 6229 since a thorough effort based on deep photometry has 
already been carried out by Borissova et al. (1999). In that paper the 
authors soundly argue from their differential age analysis that NGC 
6229 and M5 are coeval within 1 Gyr. To the best of our knowledge 
no direct estimation of the absolute age of NGC 6229 has been car¬ 
ried out. However, the age of M5 was determined by Jimenez & 
Padoan (1998) via its luminosity function and they found an age of 
10.6 zb 0.8 Gyr taking into account [a/Fe]=-i-0.4. Two more recent 
determinations of the age of M5 are by Hotter et al. (2010) using 
relative ages isochrone fitting and by VandenBerg et al. (2013) us¬ 
ing an improved calibration of the ’’vertical method” or the magni¬ 
tude difference between the turn off point and the HB; these authors 
find 12.25zb0.75 Gyr and ll.50zb0.25 Gyr respectively. Thus, in the 
CMD of Fig.j^we have overlayed the corresponding isochrones for 
12.0 Gyr from the Victoria-Regina evolutionary models (Vanden¬ 
Berg et al. 2014|^and for the metallicities [Fe/H]=—1.42 (red) and 
[Fe/H]=—1.31 (blue) and [a/Fe]=-i-0.4, shifted for the average ap¬ 
parent distance modulus /i=17.440 mag found from the RRab and 
RRc stars, and E{V - /) = 1.616 E{B - V) = 0.016. 

The isochrones for these two metalicities are very similar ex¬ 
cept perhaps at the tip of the RGB where the richer isochrone 
gets less luminous. However, if a vertical shift is applied to these 
isochrones within the uncertainty of the distance modulus, 0.01 
mag, then the two cases are indistinguishable. We can stress then 
that our CMD is consistent with the metallicity and distance de¬ 
rived in this paper and an age of 12.0zbl.0 Gyr found in the above 
papers for the coeval cluster M5. 


7 SUMMARY OF RESULTS AND CONCLUSIONS 

We have significantly updated the census of variable stars in NGC 
6229. The DIA approach to the CCD image reduction and analysis 
allowed us to recover all of the 48 variables previously known and 
to identify 25 new variables; 10 RRab, 5 RRc, 6 SR, 1 CW, 1 SX 
Phe and two new variables that we were unable to classify. 

Numerous new Blazhko variables have been identified among 
both the already known and the new RRL stars, including five RRc 
variables (V16, V17, V47, V61 and V62) (see Table|^. Also our 
data do not support the Blazhko effect seen in V5, V18, V33, V34, 
V35, V37, V46 and V48 by BCVOl based on the Dm parameter 
which is very sensitive to light curve coverage and scatter. 

Thanks to the inclusion of archive data from 1932-1935 
(Baade 1945) and from 1956-1958 (Mannino (1960) we have noted 
secular period variations in 16 RRL and have calculated the period 
change rate in 13 of them. 


^ http://www.canfar.phys.uvic.ca/vosui/#WRmodels 


The variability of V29 and V30 is confirmed and they are clas¬ 
sified here as RRab. On the other hand the variability of previously 
suspected variables V25, V26, V28 is not confirmed. The classifi¬ 
cations of V31, V32 and V39, previously believed to be EB, RRc 
and RRab respectively, are now corrected as RRab, RRab and RRc. 

The metallicity and distance of NGC 6229 were calculated 
via the light curve Eourier decomposition of RRab and RRc stars to 
find averages of [Ee/H]uvES = —1.31 zb 0.12 and d — 30.3 zb 1.4 
kpc. These results agree with previous independent determinations 
and the observed CMD matches the above results and the isochrone 
of age 12.0 zb 1.0 Gyr (Jimenez & Padoan 1998; Hotter et al. 2010; 
VandenBerg et al. 2013). 

Regarding the RRL distribution in the HB we have stressed 
once more the clear splitting between stable RRc and RRab with 
an empirical border at about (U — /)o = 0.45 — 0.46 mag, 
clearly observed in NGC 6229 and equally clear in NGC 4590 and 
NGC 5024, despite being of different Oosterhoff type. This border 
is interpreted as the red edge of the first overtone instability strip. 
Blazhko variables and double mode stars do not obey that border. 
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APPENDIX A: COMMENTS ON INDIVIDUAL STARS 

In this appendix we comment on the light curves, variable types 
and nature of some interesting or peculiar variables in TableWe 
put some emphasis on the amplitude and phase modulations of the 
Blazhko type in specific stars and note that several variables cata¬ 
logued as Blazhko variables by BCVOl are not confirmed by our 
data and that we have detected a few new ones. In their study, these 
authors combined older photographic data converted into B mag¬ 
nitudes with their own B DAOPHOT and ISIS photometry to eval¬ 
uate the amplitude variations (see their Fig. 7 for example). Then 
they calculated the deviation parameter Dm, (Jurcsik & Kovacs 
1996; Kovacs & Kanbur 1998) and for classification purposes as¬ 
sumed that RRab stars with Dm > 5 are Blazhko variables but note 
however, that examples exist of Blazhko and non-Blazhko variables 
that do not obey such criterion as Dm which is quite sensitive to 
the coverage of the light curve. In the present case, all the stars we 
catalogue as Blazhko variables are judged from our homogeneous 
data sets in the V and / filters, where the amplitude modulations 
are neatly distinguished in the light curves in Figs.[7]and[^ Hence, 
some inconsistencies in the identification of Blazhko variables be¬ 
tween the two studies are expected and we make them explicit in 
the following paragraphs. 

AO.l RR Lyrae stars 

V2, V9, V12, V43. These stars show Blazhko modulations not 
identified as such by BCVOl. Stars V2, and V9 also undergo secu¬ 
lar period variations which are reported in Table 
V3, V59. These RRL stars display very scattered light curves de¬ 
spite being located in an uncrowded region. We find no evidence of 
either star to be blended. However, we note relatively poor image 
subtractions for these stars and in this area of the field, possibly 
due to the proximity with the image edge. These are not Blazhko 
modulations. 

V5, V18, V33, V34, V35, V37, V46, V48. We find no evidence of 
amplitude modulations in these stars classified as Blazhko variables 
by BCVOl. 

V7. The star was classified as a RRab Blazhko variable by BCVOl 
based on its peculiar Fourier coefficients. Our data show very 
prominent amplitude modulations ranging about half a magnitude, 
but we note that RRab-like variations are preserved. Our data are 
insufficient to estimate the Blazhko period, but an inspection of the 
amplitudes and dates shows that the first amplitude decrease took 
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place between April 11th and July 9th 2011. A new amplitude in¬ 
crease occurred from March 5th to June 18th, 2013, i.e. in each 
case the amplitude extrema were reached within 90 and 100 days. 
V16, V17, V47, V61, V62. These RRc stars show evidence of am¬ 
plitude and phase modulations that cannot be solved by a secular 
period change or explained by double mode pulsation. We consider 
these modulations as due to the Blazhko effect. We note that star 
V47 falls much to the blue from the HB. However, it is an unre¬ 
solved pair (see Fig. and undoubtedly its colour has been pol¬ 
luted by the neighbour. 

V29. It was listed by Spassova & Borissova (1996) as a probable 
variable star and it was not included in the study of BCVOl. Our 
data plotted in Fig. |7] clearly show a RRL type light curve with 
a small amplitude modulation. Thus we confirm it to be a RRab 
variable with the Blazhko effect. For its period (0.629509d) its am¬ 
plitude in both V and I is peculiarly large, as seen in Fig.p^ 

V30. This is a faint and blended star quite close to the cluster core 
displaying a very noisy light curve. Despite this, its light curve 
shows a RRab-like shape. The star was reported as a possible vari¬ 
able by Spassova & Borissova (1996) and was not included in the 
paper by BCVOl. Hence, our light curve in Fig. is the first one 
published for this star. We confirm the variability and suggest that 
the star is likely a background RRab which would be consistent 
with its fainter and redder position on the CMD. 

V31 and V55. V31 was classified by BCVOl as an eclipsing binary 
with a period of 0.698893d. Our data do not phase well with that 
period but instead we find 0.53785Id. The light curve phased with 
this period appears like that of a RRab star (Fig. [^. V31 is very 
near a new variable star identified in this work, V55, also a RRab 
star. The periods of these two stars are similar but not equal, since 
each one does not phase the light curve of the neighbouring star. 
Also the times of maximum do not coincide. We have confirmed 
that both stars are clearly variable by blinking the difference im¬ 
ages. The amplitude modulations visible in both stars however, are 
most likely not authentic as the light curves of both stars may be 
influencing each other. Nevertheless, we do not confirm V31 as an 
eclipsing binary but instead we classify it as RRab. Both V31 and 
V55 lie well within the HB and, in the Bailey diagram (Fig. p^ , 
they also fall right on the RRab stars locus. 

V32. This star is reported as RRc in the CVSGC. We find how¬ 
ever that the period 0.603805d phases the light curve better and 
produces a RRab shape with amplitude modulations. The star’s po¬ 
sition in the Bailey diagram (Fig.|10|) and in the CMD confirm the 
above classification. 

V35, V41, V45, V46. The light curves of these RRab call attention 
for their small amplitudes which, however, are consistent with their 
long periods as can be corroborated in Fig.p^ 

V36. It is listed as a possible double mode variable in the CVSGC. 
We have found no signs of a secondary frequency in our data but 
just a period of 0.264133d. We classify this star as a RRc and it is 
labelled accordingly in the CMD and the period-amplitude plane. 
V39. This star is listed as RRab in the CVSGC. However, our data 
favour a period, light curve shape and position on the Bailey dia¬ 
gram and CMD all consistent with the star being a RRc. It must 
be stressed that the light curve in BCVOl is missing the rise to 
maximum and the maximum itself, which may have biased their 
classification. 

V48. This star has been classified by BCVOl as a RRab with a 
period of 0.516412d. Our light curve shows a poor phasing with 
that period and a glance of the light curve in BCVOl also reveals a 
poorly phased sinusoidal curve. We find that a period of 0.340270d 
phases the light curve well (Fig.[^. This period and the amplitude 


place the star well within the RRc star domain in the Bailey diagram 
(Fig.p^. Thus we have reclassified the star as RRc. 

V50. This RRab star shows clear amplitude modulations which we 
attribute to the Blazhko effect. 

V51, V58, V62, V63. By their periods, light curves morphology 
and position in the Bailey diagram of Fig.p^these are all clear RRL 
stars. However, in the CMD they share a common, peculiar posi¬ 
tion to the red and above the HB. The most extreme cases are V51 
(RRab) and V63 (RRc); both are unresolved blends (Figs. [^and[^ 
which, at least partly, explains their low amplitudes and hence their 
peculiar position in Fig.[^ V58 (RRab) and V62 (RRc) fall at the 
expected place on \hQ Ay — logP plane and are more isolated but 
perhaps not enough for their colours to be unaffected. 


AO.2 SR stars 

V64-V66, V70-V72. Their light curves as a function of HJD are 
displayed in Fig. |A1| While the variations are clear, given their 
likely long characteristic times of variation, our data are insuffi¬ 
cient to attempt a period determination. 


AO. 3 Population 11 variables 

V8, V22. The periods of these two CW stars given in the 
CVSGC are confirmed in this work (Tableand the phased light 
curves are displayed in Fig. |A2| 

V27, V73. V27 was classified as a Pop II cepheid by BCVOl 
with a period of 1.13827d; they noted however that their light curve 
was not well defined and regarded the period as uncertain. Our data 
are poorly phased with this period. Our data are plotted as a func¬ 
tion of HJD in Fig I A3 1 and they display clear variations that are 
difficult to reconcile with one single period. V73 is a newly found 
variable and, like V27, shows irregular but clear variations. Our 
data seem insufficient for a more detailed analysis of the charac¬ 
teristic time of variation in these two stars. Their positions on the 
rms and the CMD diagrams of Figs. andseem to support their 
classification as CW stars. Hence we tentatively classify them as 
such. 

V67. In the CMD this newly discovered variable is found about one 
magnitude above the HB and near V27. Its period of 1.57547d and 
light curve shown in Fig. |A2| make it a clear CW variable. It is of 
interest to comment whether or not V67 is an anomalous cepheid. 
Given the distance modulus of NGC 6229, its absolute magnitude 
is ~ —0.5. Plotted on the P — L relation for anomalous cepheids 
pulsating in the fundamental mode (see fig 6. of Pritzl et al. 2002a), 
V67 falls well below the main relation, i.e. it is about 0.5 mag 
fainter than anomalous cepheids with a similar period. Hence V67 
is not an anomalous cepheid. 


AO.4 Variable blue stragglers 

V68 and V69 show variability in our data and are contained in 
the blue stragglers region. Their light curves are shown in Fig. |A4| 
V68 is a clear SX Phe star for which only one period has been 
found, 0.0384563d. V69 on the contrary displays a small ampli¬ 
tude in V and our 1 data are too noisy to see the variation. Its pe¬ 
riod 0.27I264d is rather large for an SX Phe star thus we refrain 
to classify it as such. Alternatively, the star could be a background 
star of a different type, perhaps a RRc. 
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Figure Al. Light curves of the six new SR stars. 
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Figure A2. The three CW stars in NGC 6229 phased with their periods in Table 3. For V8 the data from Baade (1945) (orange) and Mannino (1960) (green) 
are included and show the stability of its period. In V67 small amplitude modulations are visible. 


AO. 5 Variables not confirmed 

V25, V26, V28. These stars were listed among a group of twelve 
possible variables by Spassova & Borissova (1996). The numbers 
were assigned by BCVOl although they did not include them in 
their study (as also did not include V29 and V30). We found no 
convincing signs of variability. The three stars are very close to the 
cluster core and then their light curve RMS values are large; 0.054, 
0.035 and 0.071 mag in V and 0.030, 0.024 and 0.036 mag in I 
respectively. Despite not being variable, these stars are identified in 
Figs.[^and[^and their photometry included in Table|^for the sake 
of completeness. 
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Figure A4. Two variables in the blue stragglers region. V68 is an SX Phe 
star. We have not been able to classify V69. However, given its period, V69 
could be a background RRc star, see § |A0.4| 
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Figure A3. Light variations of the possible CW stars V27 and V73. Black and red circles are the V and / magnitudes respectively. The / magnitudes have 
been arbitrarily shifted to accomodate them in the same box, so that the variations in both Liters can be appreciated. Note that V-I variations are also evident 
in both stars. 
























